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ABSTRACT 

It is now widely accepted that globular cluster red giant branch stars owe their strange abundance 
patterns to a combination of pollution from progenitor stars and in situ extra mixing. In this hybrid 
theory a first generation of stars imprint abundance patterns into the gas from which a second generation 
forms. The hybrid theory suggests that extra mixing is operating in both populations and we use the 
variation of [C/Fe] with luminosity to examine how efficient this mixing is. We investigate the observed 
red giant branches of M3, M13, M92, M15 and NGC 5466 as a means to test a theory of thermohaline 
mixing. The second parameter pair M3 and M13 are of intermediate metallicity and our models are able 
to account for the evolution of carbon along the RGB in both clusters. Although, in order to fit the 
most carbon-depleted main-sequence stars in M13 we require a model whose initial [C/Fe] abundance 
leads to a carbon abundance lower than is observed. Furthermore our results suggest that stars in M13 
formed with some primary nitrogen (higher C-l-N-l-0 than stars in M3). In the metal-poor regime only 
NGC 5466 can be tentatively explained by thermohaline mixing operating in multiple populations. We 
find thermohaline mixing unable to model the depletion of [C/Fe] with magnitude in M92 and M15. It 
appears as if extra mixing is occurring before the luminosity function bump in these clusters. To reconcile 
the data with the models would require first dredge- up to be deeper than found in extant models. 

Subject headings: (Galaxy:) globular clusters: individual (M3, M13, M15, M92, NGC 5466), stars: 
abundances, stars: evolution, stars: Population II 



1. INTRODUCTION 

It is clear from the annals of stellar astrophysics 
that low-mass stars undergo mixing as they ascend 
the red giant branch (RGB hereaf ter, for example see 
ICharbonnel fc Do Nascimentd Il998[ ). This is contrary 
to canonical models of red giants that incorporate con- 
vection as the only form of internal mass transport. 
Such models predict just one event capable of chang- 
ing surface abundances of light elements on the first 
ascent of the red giant b ranch, nam ely the so-called 
first dredge up (FDU e.g., |lbi3[i963)- 

However, a va- 
riety of observations indicate that some form of non- 
convective mixing (which is often termed "extra mixing" 
in the literature) progressively transports CN-processed 
material to the surface along the upper RGB whilst re- 
ducing the "'^^C abundance. The effects of extra mix- 
ing are observ ed in field stars (Spite ct al. 2006, 200E 
Gratton et all[2000t ICharbonnel et al.lll998: Sncdciiet^ 



burning shell progresses outwards through the star con- 
tributing processed material to the inert core. A convec- 
tive envelope which transports material to the stellar sur- 
face is separated from the shell by a radiative region. Ex- 
tra mixing requires some mechanism to transport material 
across this radiative zone to the top of the hydrogen shell 
where it encounters temperatures where burning can oc- 
cur. Processed material is cycled back to the convective 
envelope accounting for the observed compositions. Iden- 
tifying the exact physics that induces this mixing has been 
debated at length; it has proven to be a challenging task 
for stellar astrophysics. 

Phenomenological models of extra mixing can 
match the general abundance trends associated 
with further d re dge u p of CN processed m aterial 
(iPalmerini et al] 2011aB: Denissenkov &: VandenBerei 

iJli 



1986^ and globular cluster (GC) stars (iShetrond lfoOS 
Recio-Blanco fc de Lavernvl I2007D with increasing effi- 



ciency seen at lower metallicities. The mixing occurs in- 
ternally as the star ascends the giant branch therefore 
variations in abundance manifest themselves as a func- 
tion of luminosity. This is most clea rly seen in car- 
bon (iSmith fc Martel]|l2003D and lithium ()Lind et al.ll2009t 
iMucciarelh et al.ll2011l )~ 

During the RGB phase of evolution the star possesses 
a degenerate hydrogen-exhausted core. A thin hydrogen- 



l2003t Wasserburg et al.l Il995 ). Many physically based 
mechanisms have also been explored as potential can- 
didates: rotational mixing (fSwcigart & Mcnge| 119791 : 
[Chaname et al. 2005; Palacios ct al. 2006 ), niagnetic fields 
IPalmerini et al.. 2009: Nordhaus et al .200l iBusso et all 
120071: lHubb~ fc Dearbor nl ll980D . internal gravity waves 
(jDenissenkov fc Tout' '2000^ and more rece ntly thermo- 
haline mixing ( Eggleton ct al. 2006; Charb onnel fc ZahnI 
l2007at lEggleton et al.ll2008f ) and the co mbination of ther- 
moha line mixing and magnetic fields (jDenissenkov et al.l 
l2009t ). Like any new paradigm thermohaline n iixing has 
stimulated a wealth of subsequent work (see IStancliff^ 
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2010t ICantiello fc Langeil 120101: ICharbonnel fc Lagardd 
201C; lAngelou et alj 2011 ). The mechanism is initiaUy 



driven by instabihties brought about from local in- 
versions in the mean molecular weight. These inver- 
sions arise through the burning of "^He. The reaction 

■^He ('^He, 2p) He produces more particles than it con- 
sumes. If this reaction can occur in a homogenized 
region, such as is the case when the hydrogen burn- 
ing shell advances through the composition discontinu- 
ity left by first dredge- up (FDU), then its effect on the 
molecular weight profile can dominate ( for a detailed 
expla nation of the mixing process see lAngelou et al.l 
1201 ih . Understanding and modelling t his mechanism 
is an area of activ e researc h at present fWachli n et al.l 
201lHTraxler et al.i r2011: De nissenkov fc M crrvficldl l201ll: 



Charbonncl & Zahn 2007b). 

Charbonnel fc Zahn (,2007al ) have shown that the mech- 



anism can match the RGB abundances of field stars, with 
a dependen ce on metallic ty which was also demonstrated 
by Egglcto n et al.l (|2008D . The next logical step was to 



test the mechanis m against abundance measurements in 
globular clusters as lAngelou et al.l (|201lD did with stars in 
M3. M3 was chosen as a test case because it is considered 
a typical glo bular cluster. A metallicity of [Fe/H] = — 1.4 
(jSneden et al.. 2004; .Cohen fc Mel cndcz 2005) means that 
it falls near the mode in the metallicity distribution of 
halo globular clusters. It i s well studied and ex hibits sig 



nificant fC/Fe] depletion (ISuntzefj Il98lt iNorris & Smitli 
1984'; 'Smi th et al.l 119961: ISmithI 120021 : ISmithfc Martel 
12003,: Martell et al.ll2008b|) as well as redu ction in I2c/i3r 



( Pilachowski et al.|[200a 



Pavlenko et"all I2003D along the 



RGB. Here we extend this work and compare M3 to its 
second parameter partner M13 thus testing the mechanism 
in another cluster of similar metallicity. We also turn to 
the more extreme end of the GC metallicity distribution 
by investigating the clusters M92, M15 and NGC 5466 (all 
with [Fe/H]ps —2.2). This sample allows us to investigate 
the mechanism in multiple clusters at a range of metal- 
licites. The exercise is somewhat complicated by the fact 
that globular clusters possess multiple stellar populations. 

Historically globular clusters (GCs) have been utilized 
as test-beds for stella r theory. Some of th e earlies t pub- 
lished studies of GCs ()Arp et al.lll953 lA^. 1955: Sandagd 
Il953( ) coincided with the first stellar evolution calcula - 
tions (jOke fc SchwarzschildllTgsl llben fc EhrmanI [1961 . 
The initial mass spread of the stars and their presumed 
coeval nature contributed to the color magnitude di- 
agram's useful ness as a diagnostic tool for the devel- 
oping theory (ISandagd 11954 iJohnson fc Sanda"g3 Il955t 
llben fc FaulkneHl968f) . During this time GCs were consid- 
ered to be simple stellar populations; that is, the stars in 
any given cluster were assumed to be of the same age and 
composition. Hence the cluster color magnitude diagram 
was an outcome of variations in stellar mass. The first 
study to challenge the si mple stellar population hypoth- 
esis was that of Popped (| 19471 ). He measured CN band 
strength in the stars of M3 and M13 and found a CN- 
strong star amongst many CN- weak. The im portance of 
Popper's find was highlighted bv lOsbornI ()1971h who found 
CN-strong stars in M5 and MIO. These were the first clues 
that GCs possess heterogeneous C and N abundances. 

The heterogeneity of stars within GCs extends beyond 



carbon and nitrogen. Along the RGB, as well as on the 
main sequence, O, Na, Mg, and Al can disp lay star-to- 
star variations (e.g., iGratton et al.l |2004| and lKraftlll994l 
and references therein). It is now commonly argued that 
a previous generation of stars polluted the medium from 
which a distinct second population formed. The first gen- 
eration imprinted the abundance patterns from the var- 
ious stellar burning s ites, whether the y were asymptotic 
giant branch stars (iFenner et al.l 20041) . super a symptotic 
giant bra nch stars ([Ventura fc D'Antonal |2011|) , massive 
binaries fde Mink et al.' '2009D_ or massive rotating stars 
ferown fc Wallcrstcin 1993:" [Smithl [20061: iDecressin et al.l 
[20071 : ICharbonne]||2010l) . This scenario is consistent with 
the presence of C-N-O-Na abundance dispersions on the 
main sequences of GCs. 

The formation of a second generation of stars from gas 
that is rich in CN-processed material provides a paradigm 
that explains the observations of dichotomous CN band 
strengths in GC stars. The CN band strength is a use- 
ful indicator of nitrogen content of globular cluster RGB 
star atmospheres. New stars that form from the pol- 
luted gas are inherently enriched in nitroge n as well as 
^He c ompared to the primordial generation ([Norris et al.l 
Il98lh . The fact that the dichotomy in the CN band 
strength has been detected below the bump in the lu- 
minosity function (LF bump), albeit for clusters that 
are of intermediate to high metalli city, provides stron; 
evidence for the pollution scenario (iHesser fc Belli Hm 



2001 



eviaence tor tne ponution scenario messer &: tiem li^sut 
Brilev et al.lll991l:ISuntz eff fc Smith"1991':'Buonanno et al.l 



1994; Cannon et al. 1998; Cohen 1999; Brilcv fc CohcA 



iPancino et al.H2010D 



The CN bands can be used as population tracers in all 
but the most metal- p oor clusters whe re [Fe /H] < —2 (e.g., 
IShetrone et"aI1l20ld [Smolinski et ai1l2011l). since the CN 
bands are weak at low metallicity ([Cohen et a l. 2005). In 



such cases the existence of mu ltiple populat i ons ca n be 
inferred through other means. ISneden et all ([2OOO0 and 
ISneden et al.l ([1997D have shown that in metal poor clus- 
ters [Na/Fe] varies by over 1 dex and this variation is 
not correlated with RGB evolution. The temperatures 
required for processing of Na are beyond those reached 
by thermohaline mixing on the RGB. We therefore expect 
a normal (primordial) and enriched population given the 
large spread. This is a strong indication that multiple 
populations do exist in these clusters. 

Evidence of this multiple population scenario is not lim- 
ited to the observations of CN band stre ngth or varia- 
tions i n JNa/Fe] . It is di s cussed at le ngth by Catelan et al.l 
([2002f) . [Anderson et"al[ ([2OO90 and iPiottol ([20091 see ref- 
erences therein) for the main sequence (MS) and sub gi- 
ant branch. Techniques such as isochrone fitting require 
distinct *He abundances in these populations to explain 
color magnitude diagrams. Examples of such clusters 



inclu de 
l200| 



Omega Centaur i ([Piotto et al.l 
Sollima et al.l I2007D . NGC 2808 



20051: [Lee et al 



^, „ _ _ _ ([D'Antona et aL. 

2005 ULee et al.ll2005lr[piotto et al.l 120071) and 47 Tucanae 



( Anderson et al. '2009) . We note also that calculations by 
iD'Antona fc C aloi (2008) show the horizontal branches of 



various clusters can be best reproduced with the presence 
of two populations, one of which has elevated ^He abun- 
dances, which is a prediction of the pollution scenario. 
The prevailing picture is that in order to match ob- 
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servations, pollution from a primordial generation must 
be present in the cluster to provide an enriched envi- 
ronment from which a second population forms. Fur- 
thermor e extra mixing m ust b e operating in each gen- 
eration (ISuntzef f k Smi tbl 119911: iDenissenkov et al.lll998t 
IBrilev et all 11999,: ,SmitlJl2002(ra s they ascend the RGB. 
This is exactly what Smith (2002;) suggested is happening 
in M3. 'An gelou et al.l |201lt) provided further evidence for 
this hypothesis and inferred that the cluster is comprised 
of two population^^ with different ages and abundances. 
It was found that thermohaline mixing can account for 
the variation in carbon abundance of both the CN-weak 
and CN-strong stars but only for the nitrogen abundance 
of the CN-weak stars. The CN-strong stars have so much 
nitrogen to begin with that any extra mixing does not sig- 
nificantly affect the surface composition. Furthermore the 
spread in [N/Fe] is dominated by the variation that was 
present among the cluster stars before they commenced 
RGB evolution. 

In this work we examine the hybrid picture for the chem- 
ical evolution of GC giants across a range of metallicity. 
We investigate the variation of surface [C/Fe] with mag- 
nitude on the giant branch of various clusters. The ini- 
tial abundances are chosen to match the observed sub- 
giant values. We concentrate on C and N because only 
for these species do data exist over a wide range of lu- 
minosity for clusters of different metallicity (data for Li 
as a function of magnitude exists for M4 and NGC 6397, 
but the lithium values saturate over a range < 1 mag- 
nitude wher eas the [C/Fe] show va riati on over > 2 mag- 
nitudes, see iMucciarelli et al.l 1201 ll and iLind et al] l2009i 
Also measurements of lithium above the LF bump are 
scarce and many are upper limits rather than detections.). 
In the intermediate-metallicity clusters we compare M3 
to M13. We also tu rn to the metal poor regime where 
lAngelou et al.l (|2010[ ) alluded to significant problems with 
the contribution of extra mixing in M92; namely the stars 
appear to begin mixing before they reach the LF bump. 
Here we examine a further two similarly metal-poor clus- 
ters in M15 and NGC 5466 in order to help identify 
whether the behavior is due to a unique mixing history 
in M92 or is common amongst metal-poor systems. 

2. CALCULATIONS 

We use MONSTAR (the M onash version of the 
Mt. Stromlo evolution code; see ICampbell fc Lattanzid 
|2008|) to produce stellar models for the clusters. 
Our im plementation o f ther mo haline mixing follows 
that of lAngelou et al.l (1201 iD . Charbonnel fc Lagardd 
(|20T1. iSt ancliff d (l20M7^ iStancliffe et al.l mm 
and Icharbo nnel fc ZahnI (|2007al lb|) 
the formulation developed by 



th at w e 

Ubkhl (IT972h 



use 
and 



IKippenhahn. Ruschenplatt. fc Thomas! ( 1980() . where 
thermohaline mixing is modelled as a diffusive process 
with coefficient: 



(1) 



where ip — {d\n p/d\n fi)p x, S — — (91n/9/91nr)p^, 
- (ain^i/ainP), Vad (ainT/91nP)ad, V 



{d\nT/dlnP), K is the thermal diffusivity and Ct is a 
dimensionless free parameter. In fact Ct is related to the 
aspect ratio, a, of the thermohaline fingers (assumed to 
be cylindrical) by: 
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(2) 



Following lAngelou et al.l (|201ll) , we employ the empir- 
ically derived value of Ct= 1000 for the st a rs in GCs. 
It should be noted that ICantiello fc Langeil (|2010( ) and 
'Deni ssenkovl (j^OlO") prefer a lower value of Ct= 12. Their 
choice of parameter i s based on theoret i cal gr ounds and 
matches the choice of Kippenhahn et al. (119801 ) . Num eri- 
cal simulations from IDenissenkov fc Merrvfield (j201lt) and 
iTraxler et al.l (|201lh also indicate that a value of Ct= 12 
should be ad opted. This value i s unable to reproduce 
observations (jCharbonnel fc ZahnI l2007at lAngelou et al.l 
1201 ll) but must serve as a caveat to the approach taken 
in this work. We return to this point in Section [ 



2.1. Intermediate metallicity: MS and MIS 

lAngelou et al.l ([2011') showed that the chemical evolu- 
tion of the stars in M3 is consistent with the hybrid theory 
of globular cluster evolution. In intermediate metallicity 
clusters the contribution of primordial enrichment can be 
traced through measurements of the CN bands and these 
were used in M3 to infer the distinct populations. Fur- 
thermore in each population the evolution of [C/Fe] along 
the RGB was well modelled by thermohaline mixing. It is 
important to test whether this result pertains to other clus- 
ters of similar metallicity. As such we turn our attention to 
M13. M3 and M13 are a second parameter pair of metal- 
licity [Fe/H]«j - 1.4 and a ges between 11.3 a,nd 14.2 Gyr 
(IChabover et al. 1992; Jir aenez et al.l 19961: VandenBerd 
12000^ 'Sal aris fc Weisai2002t lAlves et al.ll2004D . Similari- 
ties between the two clusters exist in age, metalli city, Na-0 
anticorrelation (Sncdcn ct al. 1992; Kraf t et al. 1992), Na 
spread (|Coheniil978^.Peterson.il980; .Kraft et al.|[T992i ) and 



r-process variation (jRoederer 20111 ): it is only horizontal 
branch morphology that differs significantly. We therefore 
expect thermohaline mixing operating in multiple popula- 
tions to also model the variation of [C/Fe] (and [N/Fc]) in 
M13. 

We use data for M3 and M13 f rom a v ariety of sources. 
The M3 data are taken from ISmithI (|2002l) which is 
compiled from previous stud i es in the li terature, n amely 
ISuntzeia (119811) . ISmith et al.l (I1996D and Eii (I1999D with 
a zero point offset applied in order to homogenize the 
data. In addition to determining carbon and nitro- 
gen abundanc es as a function of absolu te magnitude, 
ISuntzefj (|1981D and iNorris fc Si^fithl (|1984 ) measured the 
CN strength for M3 giants. We do not possess measure- 
ments of CN strength for the stars in the M13 sample 
but (unlike M3) observations of [C/Fe] and [N/Fe] ex- 
ist on th e main sequen c e. For M13 we in c orporate the 
studies oflSuntzefj (119811) .ISmith et al.l (ll996D lBriley et al.l 
(200|), IBrilev et al.l (I2004D and ISmith et al.l (I2005D . Off- 
set s to allow for systematic differences have been applied 
by IBrilev et al.l (|2004 ) for all but the most recent study. 
We assumed a distance modulus of m — M = 14.43 which 



This is not a bold inference. It is expected that almost all globular clusters contain multiple populations. However this fact has not been 
observationally confirmed below the the RGB of M3. 
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is an average of various measure ments (m — M = \4,AA 
from Bucklev fc Longmord Il992l m — M = 14.33 from 
iHariid [19961, m - M ^14-47 f rom iGratton et all Il997l 
m - M = 14.48 from lReidiri997t) . The uncertainty intro- 
duced here is no larger than the measurement error of the 
observations. 

In Figure [T] we plot the determined [C/Fe] and [N/Fe] 
values against absolute visual magnitude for M13 and 
MS. Symbols are used to specify to which study the 
data corresponds; the legend can be found in the figure. 
We include two c alculations (the two solid curves) used 
to model M3 (see Angelou et al.|[2"01l[ ). The models are 
of mass M = 0.8Mq and metalhcity Z = 0.0005. In 
these calculations we run without convective overshoot, 
our mixing length parameter a is set to 1.75 and ther- 
mohaline mixing parameter Ct= 1000. In both clus- 
ters the solid black curve represents our CN-weak model 
with initial abundances Y = 0.2495, X{C) = 5.45xl0~^ 
X{N) = 1.5x10-5 and X{0) = 2.86x10^4. The initial 
carbon and nitrogen values were selected to match the 
measurements of these species in the CN-weak stars in 
MS. The solid grey curve corresponds to our CN-strong 
model with X{C) = 1.9xl0-^ X{N) = 5.5x10"^ and 
X{0) = 2.6 xlO^^ to match the CN-strong population in 
MS. As the stars are likely to have undergone CN cy- 
cling we increased our '*He to Y = 0.28 which only has a 
marginal effect on the location of the LF bump. 

Thermohaline mixing operating in two populations 
models the behavior of [C/Fe] along the RGB in both 
MS and MIS. In the case of MIS we have the added 
constraint of observations along the main-sequence. In 
this cluster the CN strong model acts as a lower enve- 
lope to the RGB observations and simultaneously accounts 
for most of the main-sequence data. However, there are 
a number of sub-giant and main-sequence stars at lower 
[C/Fe] than seen in the MS data. In order to match 
the entire main-sequence spread we have taken our CN 
strong model and reduced the initial carbon by a further 
0.4 dex. This is our dotted blue curve where we set ini- 
tial abundances to X{C) = 7.7xl0-^ X{N) = 8.8x10^^ 
and X{0) = 2.4x10^^. By design this model matches 
the most carbon depleted subgiants (except for the three 
extreme cases with [C/Fe] < — 1) but the predicted abun- 
dances along the RGB are systematically lower than the 
current observations. It appears as if these most-carbon- 
depleted subgiants lack counterparts on the RGB. This 
could be the role of small statistics or the effect of system- 
atic offsets from differ ent data sou r ces; r ecall we do not 
apply an offset for the ISmith et al.l ((2005) study. We are 
unable to account for the two most carbon-poor stars in 
MIS ( [C/Fe] fti -1-5)- T hese stars are at the faint end 
of the 'Bril ev et al.l (|2004f ) catalogue and hence have large 
associated errors (Cohen 2011, private communication). 

Compared to MS, there appears to be a larger star- 
to-star variation in nitrogen in MIS; a significantly N- 
enhanced model is required to reproduce the upper en- 
velope of the data. In the context of the current paradigm 
the enrichment is proposed to be due to ON cycling in 
the polluting generation of stars. To investigate this we 
constructed an ON-cycled model (dashed green line) with 
half the oxygen converted to N, giving X{C) = 1.9x10"^, 
X{N) = 1.85x10-4 and X{0) = l.SxlO-^. The figure 



shows that the change in [N/Fe] produced by thermoha- 
line mixing is negligible. The initial nitrogen in these stars 
is sufficiently high that any nitrogen brought to the sur- 
face via thermohaline mixing makes little difference to the 
envelope composition. 

We note from Figure [T] that even the ON cycled model 
(green dashed line) fails to account for the N-rich stars 
which show 1 dex more N than this model. Even process- 
ing all the O into N would not produce a sufficiently high 
N abundance. We must consider the possibility, therefore 
that C-hN-fO is higher in MIS than in MS. This is a 
testable prediction. 

As we have stated there are a la rger number of stars in 
MIS that are enriched in nitrogen; ISmolinski et al.l (|201lt) 
find the cluster contains four times as many CN-strong 
stars than CN-weak ones. In terms of the hybrid the- 
ory this is interpreted as the cluster containing a ma- 
jority of second gener ation stars. In ISmithI (|2002D and 
ISmolinski et al.l (|201lD the observations suggest that the 
MS stars are evenly split amongst the CN-weak and CN- 
strong populations. The inferred ratios of the stellar popu- 
lations within the clusters, determined from measurements 
of the CN bands a g ree w ith the theoretical predictions of 
iD'Antona fc Caloil (|2008( ). They have independently de- 
termined that more than 70% of stars in MIS and 50% in 
MS are required to be enriched in "^He (a product of CN 
cycling along with N enrichment) to reproduce the respec- 
tive horizontal branch morphologies. 

In MS there is a distinct lack of CN-strong star s at low 
luminosity. As discussed in Angel ou et al.l (|201lt) this is 
an artifact of the original [Suntzefj ( 198lh studv in which 
the lower lumin osity stars that were o bserved happened 
to be CN-weak. iNorris fc SmithI (|1984D showed that CN- 
strong giants do exist in MS at lum inosities correspond- 
ing to_JhB_^int_Jimh_of the |Su^z^^ (1981) survey, as 
did ISmolinski et al.l (|201lD . The latter show the CN band 
strength remains clearly dichotomous in MS even at low 
luminosity. We assume this will be refiected in the [N/Fe] 
abundances and therefore [C/Fe]. So whilst we possess 
measurements of the CN bands for our MS sample and a 
few observations of [C/Fe] at low luminosity, conversely 
in MIS we possess many low-luminosity observations of 
[C/Fe] but no measurements of the CN bands. Once again 
the work of ISmolinski et all (2011) provides some insight. 
They demonstrate an obvious dichotomy in the CN bands 
along the RGB of MIS. The dichotomy is not as clear be- 
low the sub giant branch as t he temperatures are too high 
to allow molecule formation. ISmolinski et al.l (|201ll ) also 
show a greater spread in the CN band values especially in 
the CN strong stars; this is refiected in the large spread in 
[N/Fe] observed in the cluster. 

2.2. Low Metallicity 

2.2.1. NGC 5466 

Carbon d epletion along the RGB o f NGC 5466 has been 
observed by 'Buonann o et all (|1985D . iFekadu et all (|2007t) 
and Shctronc et al. ( 2 0101) . In this study we use the data 
from Shct rone et al.l (120101 ) as the observations are from a 
single instrument, they cover a large luminosity range and 
for each star [C/Fe] and the CN band strength were de- 
termined. As was the case for the intermediate-metallicity 
clusters, we plot [C/Fe] as a function of absolute visual 
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magnitude. In Figures ^jji. and ^jp the NGC 5466 data 
are plotted according to CN band strength; open cir- 
cles denote CN-weak stars, filled circles denote CN-strong 
starf0 whilst triangles represent stars where the CN band 
strength is unknown. The black solid curve in Figure [2^ 
is our thermohaline mixing model with initial abundances 
of r = 0.25, X{C) = 1.6 X 10-^ X{N) = 5.03 x lO'^ and 
X{0) = 4.57 X 10-5. The model is of mass M = O.SMq 
and metallicity Z = 0.0001 corresponding to the ob- 
served [Fe/Hl= —2.2 and age between 1 1 and 16 Gyr 
(|di Cecco et al.ll2010l:lGrundahl et al.ll2000D . We run with- 
out convective overshoot, set the mixing length parameter 
to a = 1.75 and set the thermohaline mixing parameter to 
Ct— 1000. Highlighted also are the locations of the impor- 
tant mixing events along the RGB. Common to all panels 
are the solid vertical line and dotted line. The solid verti- 
cal line represents the end of first dredge-up as calculated 
in the models; this occurs at My « 0.56. The dotted line 
co rresponds to the LF bump at [Fe/H]= —2.2 according 
to iMartell et al.l (|2008b[ ) based on a metallicity-LF bump 
relation; this occurs at My ss —0.55. For each individual 
cluster we also include the dashed line which represents 
a photometrically derived value of the LF bump magni- 
tude. In NGC 5466 this was determined by iFekadu et al.l 
(|2007t ) and occurs at My w -0.2. The end of first dredge- 
up marks the end of any surface composition changes ex- 
pected by canonical theory. For stars of this mass and 
metallicity there is no visible change to the envelope com- 
position from first dredge-up: the change in surface mass- 
fraction of carbon is 6 x 10~^ in our models. The position 
of the LF bump represents the earliest point at which ther- 
mohaline mixing can begin to operate. It corresponds to 
the point when the hydrogen burning shell meets the com- 
position discontinuity left behind by first dredge-up. 

Figure ^ suggests that our solar scaled model (black 
curve) is a good fit to observations of [C/Fe] in NGC 
5466. Surface depletion in the cl uster appears to begin 
afte r the LF bumps determ ined bv lFusi Pecci et all ()1990[ ) 
and IMartell ct al. (2008b). Wc find that the bump deter- 
mined by MartcU ct al. ( 2008b) is a more c o nvincing fit to 
the m odels than that of lFekadu et all (|2007D . iFekadu et al.l 
(|2007t ) also inferred the LF bump magnitude by isochrone 
fitting and o btained a value 0.3 ma gnitudes brighter agree- 
ing with the IMartell et"all (|2008bD prediction. The deter- 
mination of the LF bump at low metallicity is difficult. 
We return to this point in Section IXTl 

Although a single CN population undergoing thermo- 
haline mixing is a satisfactory fit to the data, in Figure 
|2}3 we provide an equally plausible fit drawn by eye. The 
black solid line (which applies to Figures HJi & HF) 
tries to take into account our belief that carbon depletion 
should not occur until after first dredge-up but does not 
take into consideration any preconceptions of the mixing 
profile. If we assum e the bump luminosity calculated by 
IFekadu et al.l (|2007l ) is correct then depletion begins after 
the LF bump, as expecte d. If, however , we use the bump 
luminosity calculated by IMartell et al.l (|2008b( ) then the 
mixing may be argued to begin before the bump is reached. 

NGC 5466 requires almost no spread in the initial com- 
position to account for th e spread of carbon. In a ddition to 
measurements of [C/Fe], IShetrone et al.l (|2010[ ) analyzed 



the CN bands of the stars in their sample. The results 
hinted that a CN dichotomy, which is typical in globular 
cluster red giants, may be present but as we can see from 
Figured^ the stars do not s eparate as clear ly as those in 
the more metal rich clusters. IShetrone et al.l |2010) used a 
line ar least-squares fit in the plane of the S(3839) CN in- 
dex (jNorris et al.lll98ll ) versus absolute V magnitude to di- 
vide relatively CN-strong stars (those above the line) from 
CN-weak stars (those below the line). As i n other stud- 
ies of low- metallicity globular clusters (e.g.. IMartell et all 
l2008al) . the mean separation in CN band strength between 
the two groups is not large relative to the scatter within 
each group. Whilst the data suggest that there is a CN 
bimodality in NGC 5466, its low metallicity reduces the 
effectiveness of CN band strength as a marker of multiple 
populations in this cluster. 

2.2.2. M92 

The M92 data in Figure s |2t and |2Ji are those adopted 
by ISmith fc MarteO (|2003D and comprise data from var- 
ious sources to which offsets have been applied in or- 
der to remove systematic differences in abundance scales. 
These origin al sources are studies byiC arbon ct al. (1 9821 ). 
iLanger et all (^198^ and iBellman et al.l ()200lf ) and we 
highlight the target stars in each catalogue (see the leg- 
end in the figure). Even with typical errors of ±0.2 dex 
there is a much larger spread in [C /Fe] at any given mag- 
nitude than was seen in NGC 5466. In these M92 panels 
the dotted li ne is the photome trically derived LF bump 
identified bv iNataf et al.l (|201ll ). These authors have un- 
dertaken high resolution HST observations to determine 
the LF bump in many globular clusters. Their mea- 
sured value of the magnitude of the bump My w 0.016. 
(Vfc„,„p = 14.66 6 ± 0.013, (m - M)V = 14.65) differs 
from the older iFusi Pecci et al.l (jl990[ ) value of My ^ 
—0.39. We do not possess measurements of CN band 
strength for the M92 s t ars pl otted in Figures [21: and | 2li but 
iNorris k Pilachowskil (|1985D and ISmolinski et al.l J201lt) 
have provided evidence that a bimodality may exist. As is 
the case for NGC 5466 the separation in the populations is 
marginal, a property that appears common at low metal- 
licity. Even without the C N bands the spread in [Na/Fe] 
(jSneden et al.ll2000l [T997[) . and the spread in [C/Fe] at 
low luminosities is evidence for mixing between the two 
populations. 

The deplet ion of [C/Fe] in M92 has previously been in- 
vestigated bv lDenissenkov fc VandenBerd (|2003r ) but over 
a smaller luminosity range. Their (canonical) extra mix- 
ing formulation requires two models, one depleted by 0.1 
dex and the ot her by 0.5 dex i n [C/ Fe] to match the abun- 
dance spread. lAngelou et al.l (j2010D considered the deple- 
tion of [C/Fe] with a different implementation of thermo- 
haline mixing to that used here. They highlighted that 
M92 appeared to show surface depletion of carbon before 
FDU, a property in disagreement with not just thermoha- 
line mixing but also canonical stellar theory. This strange 
behavio r has been discussed by (at least) Mart ell et al.l 
(l2008bD.lG ratton et all (|2004D IBellman et all (|200lD . a"^ 
ILanger et al. (198^. 

This apparent pre-FDU depletion is highlighted in Fig- 
ure |2ji. Here the black solid line is an eye fit based on 



The solid star is CN strong but it is most likely a CH star and hence a binary. 
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the assumption that carbon depletion should not occur 
until after FDU while the grey curve is a best fit by eye 
to the data without any preconceived luminosity for the 
onset of carbon depletion. If we follow the latter fit, then 
by the time the end of FDU is reached the stars have de- 
pleted almost 0.5 dex in [C/Fe]. This is unusual because, 
as we have stated, at this mass and metallicity we ex- 
pect changes due to FDU to be of the order AXsurfacc(C) 
= 6 X 10~®. Between the end of first dredge-up {My 
« -fO.5) an d the determined loc ation of the LF bump 
according to lMartell et al.l (|2008bD (My w -0.5; the ear- 
liest point at which thermohaline mixing is expected to 
begin) the stars in M92 are depleted by roughly a fur- 
ther 0.3 dex in [ C/Fej. Depletion i s present but not as 
pronounced if the iNataf etall (|2011h bump magnitude is 
used. It is possible that these are the observable effects of 
other extra-mixing mechanisms and this may also explain 
the carbon depletion prior to the end of first dredge-up. 
If true, this would have significant implications for stellar 
theory. First of all such deplet i on was not seen in o ther 
clusters (jSmith fc Martelll 120031: IShetrone et al.ll2010[ ). In 
most scenarios, extra mixing is inhibited until the star 
reaches the LF bump and the advancing H-shell removes 
the molecular weight discontinuity left behind by the re- 
ceding convective envelope. In the case of M92 some stars 
on the giant branch have already depleted their [C/Fe] by 
about 0.8 dex when the models reach this stage. If we have 
to postulate that some form of mixing begins sufficiently 
early to produce this depletion, then the mixing must nec- 
essarily remove the abundance discontinuity that is itself 
responsible for the LF bump. If the apparent behavior 
is real, we believe this is a significant problem for stellar 
astrophysics. 

If we assume that M92 displays an initial spread of 
[C/Fe]=0.5 dex then can we remove the issue of mixing 
before the end of FDU? The black curve in Figure [2t is 
the same solar scaled model applied to match the stars in 
NGC 5466. As we now require a greater spread in car- 
bon to match this cluster we include a second model (grey 
solid curve) run with the same physical parameters but 
a composition that is reduced in C and enhanced in N. 
For this model we set Y = 0.25, X{C) = 5.46 x 10-^ 
X(N) = 1.60 X IQ-^ and X{0) = 4.57 x lO^^ giving 
[C/Fe]=-0.5. We have not changed Y in the model. The 
main effect of changing Y is to alter the the location of 
the LF bump and as we have seen in M3 and M13 this 
is negligible. These two thermohaline models with large 
initial carbon spread remove the apparent pre FDU deple- 
tion. Assuming such a spread in [C/Fe] raises the issue 
that there are no stars observed with [C/Fe] =-0.5 below 
a magnitude of My ~ 2. Only four stars have been ob- 
served at such magnitudes and all have [C/Fe] « 0. A 
large spread in [C/Fe] is observed in M15 and M13 below 
the sub giant branch. We expect this is the case in M92 
but require a targeted study to confirm this assumption. 

A large spread in [C/Fe] does not solve all the prob- 
lems associated with M92. Mixing still appears to be oc- 
curring before the LF bump, even for the model with an 
initial [C/Fe] =-0.5. There is approximately 1 magnitude 
between the location where the carbon abundances turn 
down and where the models suggest depletion should be- 
gin. Of greater concern is the inability of the models to 



match the upper RGB abundances of stars in M92. The 
initial abundances were selected to cover the spread in 
[C/Fe] before the onset of extra mixing. No stars that 
have depleted [C/Fe] (magnitude brighter than My = 0) 
fall along the evolutionary path predicted by the solar 
scaled model (black curve). If such stars do exist and 
are governed by thermohaline mixing, we expect them to 
be bright members in the cluster (ascending the RGB). 
It would be perverse if by chance they were not selected 
for any of the three high resolution spectrographic studies 
(unless of course, they are all in the centre of the clus- 
ter). The initially depleted model (grey curve) provides 
an upper envelope to the mixed RGB data, however the 
composition was selected to form a lower limit of the ini- 
tial M92 abundance spread. It is sobering to see how many 
data points fall below this curve. Although the two models 
with an initial spread of 0.5 dex in [C/Fe] can encompass 
fainter stars (My > 0.5 there are very few stars between 
the curves once carbon depletion has begun. Thermoha- 
line does not deplete carbon to the levels seen in M92: the 
models are unable to simultaneously account for both the 
sub-giant branch and RGB data. 

The M92 sample is a combination of three different stud- 
ies and although attempts have been made to homogenize 
the data, systematic differences in the studies still may 
lead to offsets of ~ 0.3 dex in the [C/Fe] abundances. At 
such low metallicity one has to wonder if inhomogeneity 
in the data may be affecting our interpretation of the mix- 
ing history. Extensive measurements of [C/Fe] at various 
magnitudes in this cluster would provide a robust test for 
stellar evolution models. 

2.2.3. M15 

[C/Fe] as a function of magnitude is plotted for M15 
in Figures [2J3 an d [2f. The da ta are a combination of 
those from iTrefzger et al.l (|1983l crosses) and .Cohen et all 
(|2005l diamonds) and plotted as given in these sources. 
No offsets have been applied to correct for systematic er- 
rors so the data are marked according to their source. 
As is the case with M92 we do not p ossess CN band 
streng th measurements for our sample but lSmolinski et al.l 
(|2011[ ) show there would be little separation between 
the mean values of the CN-strong and CN-weak stars, a 
common aspect of the three metal poor clusters studied 
here. The same is also true for (at leas t) NGC 5053 a 
cluster with metallicity [Fe/H] = —2.3 (jSmolinski et al.l 
I20TI . In Fi gures and |2f the dotted vertical line rep- 
res ents the magnitude of the M15 LF bump determined 
bv INataf etUT (|2011h which occurs at My « -0.075 
{Vbump = 15.315 ±0.021, (m - M)y = 15.39). This is 
in comparison to IZoccah et all (|1999D who also used HST 
data to determine a bump magnitude of My ~ —0.42. 

In Figure I2J3 we plot the two thermohaline mixing mod- 
els applied in M92. In both clusters thermohaline mixing 
does not deplete carbon to the degree seen in the obser- 
vations. Like M92 there is evidence that Ml 5 is mixing 
between the end of FDU and the LF bump. This is high- 
lighted by Figure |2f where we plot our lines drawn by eye 
rather than the stellar models. This effect would be ex- 
acerbated h ad we adopted the l ocation of the LF bump 
according to IZoccali et al.l ()1999l ) . It is possible this is an 
artefact of combining inhomogeneous data. There is some 
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luminosity overlap between the two studies and it is only 
when they are combined that it appears as if pre bump 
depletion is occurring. This same argument can be used 
to explain any pre-FDU mixing alluded to by the grey line 
of best fit in Figure A homogeneous set of observa- 
tions over this luminosity range (as exists for NGC 5466) 
is required. 

While similarities exist there are also slight differences 
between M15 and M92. A large spread of [C/Fe] is present 
in M15, which extends below the sub giant branch. This is 
justification for the inclusion of our initial [C/Fe] = —0.5 
model (grey curve in Figures [2t and [5^) . In M92 it is the 
dearth of observations at low luminosity that raises many 
questions about the behavior of the cluster and its ini- 
tial abundance spread. We note there seems to be a lack 
of stars with high [C/Fe] just before FDU in M15. We 
assume these stars exist and the current void is a result 
of combining two studies that focus on stars at different 
stages of evolution. Perhaps the most curious aspect of 
M15 is the behavior of the stars at My « — 1. The stars 
appear to suddenly increase their carbon abundance. We 
await confirmation of this behavior before speculating on 
its cause as it is inconsistent with our current understand- 
in g of stellar n ucleos ynthesis . 

Lind et all (|2009D and iMucciarelli et al.l ()201l[ ) use 



lithium to trace extra mixing along the giant branches of 
NGC 6397 and M4 respectively. Such studies provide a 
complimentary diagnostic to the variation of C and N with 
luminosity. Lithium is very sensitive to nuclear burning: it 
is destroyed at low temperatures. Any transport of mate- 
rial into warmer regions will be reflected through a deple- 
tion in the surface lithium abundance. In the case of M15 
and M92 this would hopefully help identify the magnitude 
at which extra mixing begins in these clusters. Given the 
behavior of M15 and M92 the point at which surface de- 
pletion begins (after first dredge-up) must be considered 
distinct from the LF bump (even though there are good 
theoretical reasons why the two should coincide as they 
seem to at higher metallicity) . For our purposes we can 
only rely on high level photometry (and statistical analy- 
sis) of the clusters to distinguish the location of the bump 
and, as we discuss in section [5t1 this can be difficult. 

3. DISCUSSION 

3.1. Comparison across metallicity 

In our two intermediate-metallicity clusters the spread 
in [C/Fe] is explained by distinct populations which clearly 
separate according to CN band strengt h. This is not 
the case in low-metal licity clusters (Shc trone et al.l [20101 
ISmolinski et al.ll201l[ ). Any CN bimodality is marginal 
at best. The difference between the mean value of a CN 
index such as S(3839) between CN-strong and CN-weak 
stars is so small that their distribution can be interpreted 
as a single CN population with scatter. There is still clear 
evidence for multiple populations in the metal poor clus- 
ters from the large spreads in O, Na, Al abundance within 
them. We therefore expect that a polluting generation will 
enrich the cluster with CN processed material and '^He. 
This should be reflected in the carbon and nitrogen abun- 
dances. In M15 there is a la rge spread in both [C/Fe] (w 
1.5 dex) and [N/Fe] (« 3 dex,|Cohen et al."2005') but only a 
marginal change in the CN bands (,Smolinski et al..,2011i) . 



In these low metallicity clusters the CN band strength is 
not necessarily representative of the nitrogen abundance 
in the stars. 

According to our stellar models with thermohaline mix- 
ing, the clusters studied here are close enough in metallic- 
ity that we should expect a similar degree of carbon de- 
pletion. In the intermediate metallicity clusters (M3 and 
M13), we are able to match the depletion of [C/Fe] along 
the RGB using two models of different initial abundances. 
In the metal poor systems thermohaline mixing can ac- 
count only for the evolution of [C/Fe] in NGC 5466. In 
M15 and M92 not only do the clusters appear to have de- 
pleted far more carbon than we predict but this depletion 
has begun before the LF bump. Such behavior is not only 
inconsistent with thermohaline mixing but also standard 
stellar evolution. 

It is worth noting that the LF bump in extremely metal- 
poor clusters (e.g., M92, M15 and NGC 5466) is not as 
clearly visible as it is in more metal-rich clusters. There 
are two theoretical reasons why this is so. Firstly the depth 
of first dredge-up is a function of metallicity. Metal-rich 
stars have deeper convective envelopes than metal poor 
stars and therefore a greater discontinuity in the hydrogen 
profile that develops after first dredge-up. The burning 
shell in metal-poor stars will encounter a smaller hydro- 
gen difference and spend less time readjusting the stellar 
structure. Because less time is spent at the magnitude of 
the bump the likelihood of observing stars at this location 
is reduced. Secondly, evolution on the RGB speeds up as 
stars move towards brighter luminosities, so the expected 
number of stars scales inversely with the luminosity in a 
volume limited sample (but not in a magnitude limited 
sample). Hence the higher the luminosity at which the 
bump occurs (i.e., the lower the metallicity), the lower the 
overall number of stars one should expect to observe and 
t he harder it is to i dentif y the bump. 

iFusi Pecci et al.l (|1990[ ) identified an LF bump in M92 
by co-adding data for three very similar clusters. Their 
determin ed valu e of Afy « —0.39 differs slightly from that 
of Mar tell et al.l ( pOOSb ) who find My « -0.55 based on 
the use of multiple clusters to determine a L F bump- 
metallicity relation. Work by iPaust et"aLl ([2003) provides 
little evidence for a b ump in the observe d LF of M92 how- 
ever recent work by iNataf et al.l (|2011l) using HST data 
suggests that i t is present at a mag nitude of My « 0.016. 
For NGC 5466 lFekadu et al.l (|2007[ ) identified a bump us- 
ing statistical arguments. The bump in M15 was more 
easily i dentifiable through the use of HST data a s was 
done bv lNataf etall (|201lD and lZoccah et al.l (| 19991) . 

One possible cause for altering the location of the bump 
(and hence the onset of thermohaline mix ing) in our mod - 
els is to decrease the ''He abundance fSweigart| 119781 ). 
However, extrapolating from the Sweigart (1978,) results 
to match the observations would require an initial '*He 
well below the Big Bang Nucleosynthesis value. Alterna- 
tively, if first dredge-up penetrates deeper than the mod- 
els predict, then the hydrogen burning shell will encounter 
the homogenized region at lower luminosities. In Figure 
[3| we plot our metal-poor model with solar scaled abun- 
dances, in which the convective envelope extended down 
to a mass of M = 0.368 Mq. We have included a model 
(grey line) with the same initial abundance but artificially 
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homogenized (mixed from the surface) down to a mass of 
M = 0.320 Mq at the time of FDU. The results show that 
extending the depth of first dredge-up by about 15% in 
mass translates to the mixing beginning about one mag- 
nitude fainter. Note that the artificial model undergoes a 
period of adjustment as it returns to the structure of the 
standard metal-poor model. During this brief phase the 
burning shell is closer to the homogenized region than the 
structure would dictate if deeper dredge-up occurred nor- 
mally. As a consequence the temperatures are such that 
■^He is easily destroyed but carbon deplet ion is stil l inefh- 
cient; this can be seen in Figure |31 Ang elou et al.l (|201lh 
call this situation a "burning limited" regime. The mate- 
rial is efficiently transported to the desired regions but the 
processing is inefficient due to the burning conditions. The 
artificial model therefore has less "^He to drive the mixing 
near the tip of the RGB when the co nditions are mo r e con - 
ducive to carbon depletion, a regime [Angelou et al.l (j201lh 
call "transport limited". In this scenario the inefficient 
transport of material to the burning regions is the limit- 
ing factor in the processing. This however is a side issue 
stemming from the artificial model; the important point is 
that one way for the models to match the observations is 
to make first dredge-up occur deeper in the low-metallicity 
stars but only low-metallicity stars in M92 and M15. 

3.2. The inconsistency at low metallicity 

We are left wondering why there is such an inconsis- 
tent picture of globular clusters at low metallicity. M92 
and M15 possess large spreads in [C/Fe] whilst NGC 5466 
shows a much narrower distribution. We find that in the 
cluster where we have uniform observations thermohaline 
mixing appears to model the depletion of carbon well. We 
have already discussed the possibility that combining stud- 
ies in M15 may be effecting our interpretation of the data. 
Although data for M92 has been homogenized it is still the 
amalgamation of three different studies. Compiling data 
from various sources may lead to systematic errors. Sys- 
tematic offsets of 0.3 dex in [C/Fe] are possible and could 
be the cause of the apparent contradiction. Because we 
have two clusters at the same metallicity in M92 and M15 
displaying similar behavior one has to wonder how likely it 
is that systematic offsets conspire in such a way to appear 
to give the same effect. 

We note that the clusters showing the largest discrep- 
ancy with standard models, are also the most massive. 
Having measured th e integrated V magn it udes o f the three 
metal-poor clusters iGnedin fc Ostrikeil ()1997l) find that 
NGC 5466 (w 1.33 x 10^^ Mq) is a lower mass system 
than M92 (?s 3.64 x 10^ Mg) and M15 (« 9.84 x 10^ Mq). 
This raises the question whether there is a correlation be- 
tween cluster mass and the extent of primordial [C/Fe] 
spread at low metallicitites. The clusters M15 and M92 
may have been able to sustain a greater degree of primor- 
dial inhomogeneous carbon enrichment than NGC 5466, 
with a resulting greater dispersion in [C /Fe] at all magni- 
tudes on the main sequence and red giant branch than in 
NGC 5466. The greater range in [C/Fe] along the RGB 
of M15 and M92 can contribute to a greater observational 
uncertainty in identifying the magnitude at which carbon 
depletions produced by extra mixing set in. By contrast, 
a smaller primordial carbon spread in NGC 5466 produces 



a more tightly defined locus of [C/Fe] versus My on the 
RGB of that cluster, making it a preferable test case for 
mixing studies. When comparing the r-process elements 
both M92 and M15 show similar star-to -star variations 
(|Roedereill201lt iRoederer fc Snedenll201lD fa lthough this 
has recently been questioned bv lCohenI (l20Tlh ') which may 
add additional support to the Na and O abundance evi- 
dence that both clusters have sustained heterogeneous en- 
richment across a range of chemical elements. It would be 
interesting to determine the O, Na, and r-process element 
patterns in NGC 5466 to see if they are more homogeneous 
than in M92 and M15, as a mass-dependent primordial en- 
richment scenario would anticipate. 

3.3. The role of thermohaline mixing 

It is of course possible that thermohaline mixing does 
not govern the surface composition of low mass giants. 
It is not the only process by which mixing may occur in 
the radiative zone. Mechanisms such as those listed in 
Section [T] may also be involved. The way in which these 
mechanisms interact is uncertain and as they are inher- 
ently three-dimensional our understanding of their behav- 
ior will improve w ith the study of these processes in hydro- 
dynamical codes (iDearborn et al. 200lUBazan et aLl2001 



Turcotte et all 120021: lEggleton et all 120021: iBazan et all 
20031) . Still, the current generation of ID codes are 



the precursors to more sophisticated modelling and much 
insight can be gained through them. One-dimensional 
spherically-symmetric models with multiple extra- r nixin; 
mec hanisms include th ose of ' Cantiello fc Langeil ()2010[ 
and lCharbonnel fc Lagar de ( 20l"of), both of whom demon- 
strate that the thermohaline mixing diffusion coefficient 
is larger than the radial component of rotational mixing. 
The former also show the thermohaline coefficient to be 
larger than that of magnetic buoyancy. These codes do not 
explicitly treat the interaction of the mixing mechanisms 
and how any given instability may react due to the pres- 
ence of other processe s; rather, they simply ad d the diffu- 
sion co-e fficients. Cha rbonnel fc ZahnI (l2007b[ ) have inves- 
tigated an instance of multi-process interaction through 
linear analysis and argue that magnetic fields could serve 
t o inhibit the effects of thermoh alin e mixing. 

ICharbonnel fc Lagardd (|2010D and lLagarde et al.l (|20lH) 
have produced detailed rotational-thermohaline mixing 
models for stars of various mass and metallicity. They 
find rotation leads to a deeper penetration of the convec- 
tive envelope than their non-rotating models. This allows 
rotational models to begin thermohaline mixing at fainter 
magnitudes than their static models; this may go some 
way to reconciling the magnitude at which stars in M92 
and M15 begin mixing. It should be noted that low-mass, 
low-metallicity rotati onal models have been produce d by 
iPalacios et al.l ()200 6') and Denisse nkov et all ()2006[ ) and 
in the former, the depth of first dredge-up depends on the 
rotational model used. A window of 0.13 magnitudes is 
possible for the location of the bump depending on the 
treatment of rotational physics. That is, different imple- 
mentations of rotation can move the location of the LF 
bump to higher or lower magnitudes. Different models 
of rotational mixing are at play here, in much the same 
way that different codes produce different third dredge- 
up results based on the treatment of convective stability 
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(jHerwigj [20051 ). However stochastic variations such as ro- 
tation, may lea d to deeper firs t dredge-up in the stars of 
M15 and M92. iPetersonI ()1983[) have shown that stars on 
the horizontal branch of M13 rotate twice as fast as those 
on the horizontal branch of M3. Although we find no ma- 
jor discrepancies between the locati on of the LF bum p in 
these two clusters (0.07 magnitudes iNataf et al.ll20TTI ). in 
the case of M15 and M92 the effect of rotation on the depth 
of the FDU may be more pronounced. Why the stars in 
M13 would rotate faster than the stars in M3 would also 
require an explanation. 

Whilst it is generally agreed that the burning of 
^He in a homogenized zone causes an instability, the 
efficiency of the resultant mixing has been the cause 
for debate. The exact value of Ct to adopt re- 
mains contentious. Recent multi-dimensional mod- 
els of thermohaline mixing by [penissenkoy ( gOlO) , 
iDenissenkov fc Merrvfieldl (l201lh andlTraxler et~all (|201l[ ) 



support the view of ICantiello fc Langeil (|2010D and the 
suggestion by iKippenhahn et al.l ( T980f ) that the "aspect 
ratio", a, of the rising element should take a value of 
a ss 1 (CfW 12). La borato r y exp eriments by Stommel and 
Fall er published in ISterij (ll96C)D have lead lUhichI (|1972| ) 
and ICharbonnel fc ZahnI (|20073) to prefer a value of a 
somewhere closer to a w 6 (CtW 1000). Indeed there 
is great empirical evidence to prefer this value as it ap- 
pears to match observations of stars across a great range 
of mass and metallicity. The mechanism is an elegant 
means of matching C and N abundances in various stars 
as well as ensuring measurements of '^He in HII regions 
are consistent with predi ctions from stella r models and 
Big Bang nucleosynthesis (i Rood et al. 1984 : HoganI 
Charbonncl 1995; Charbonnel fc Do Nascimentol 



1995 



19981 



Tosiiil998i : iPalla et al.ii20Q0t iRomano et al.ll2003[ ). Canoni 



cal models predict that low-mass, main-sequence stars are 
net produ cers of '^He which is returned to the ISM through 
mass loss. iHata et al] (|1995D have shown that about 90% 
of the ^He produced on the main sequence must be de- 
stroyed to reconcile the two fields. Our preferred value of 
Ct— 1000 matches the carbon and nitrogen in clusters and 
also destroys over 90% of the '^He in the star before the tip 
of the RGB. The lower value of Ct= 12 once again creates 
a discrepancy between measurements of ^He in HII regions 
as only sa 20% of the initial '^He is destroyed meaning low- 
mass stars are once again produc ers of '^He. 

It is the aim of this work and in lAngelou et al.l (|2011f) to 
rely on empirical evidence to determine the best value of Ct 
to use in the ID codes. In the models presented here Ct has 
been fixed at a value of 1000. This is consistent with previ- 
ous studies where the same fact or has reproduced ob serva- 
tions of globular cluster stars lAngelou et al.l 1201 iD . field 



stars ([Charbonnel fc Lagardell2010D and where the mixing 
has been used to explain the dichotomy between extremely 
metal poor stars and carbon enhanced metal poor stars 
(jStancliffe et al.ll2009D falthough the latter po i nt is in dis- 
agreement wit h Denissenkov fc Pinsonneau a ((2008lbh l 
lAngelou et all (|20lH ) demonstrate that once Ct exceeds 
1000, the processing enters into the burning limited regime 
whereby increasing the diffusion coefficient has little ef- 
fect on the final surface abundances. They also show val- 
ues of Ct < 600 simply lead to less mixing than is re- 
quired to match the observations. Although we are mod- 



elling the mixing as diffusive it is in fact an advective 
process. We have implemented a linear model of ther- 
mohaline mixing which operates via a diffusion equation. 
Recent 3D numerical simulation s of thermohaline m ixing 
(IDenissenkov fc Merrvfieldl 120111 : iTraxler et a"ni201lD find 
blob-like structures, which are identified with the aspect 
ratio of the idealized "fingers" in the ID derivation of 
lUhichI 11972). It is from this identification that the pref- 
erence for low values of Ct in the numerical simulations 
is based. Numerical simulations are a very powerful tool, 
but are subject to numerous caveats - resolution (time and 
space), viscosity (numerical and molecular), and extrapo- 
lation to stellar conditions (e.g., Prandtl numbers that are 
« 10^ too large). In view of these uncertainties we feel it 
premature to draw definite conclusions from the numerical 
simulations. Nonetheless, if the low values of Ct coming 
from the simulations prove to hold then it would result 
in a conclusion that either the diffusive approximation to 
thermohaline mixing used here is not applicable or that 
thermohaline mixing is too weak to produce the observed 
carbon depletions of globular cluster red giants. 

4. CONCLUSIONS 

We have modelled the depletion of carbon via thermo- 
haline mixing on the RGB of the intermediate metallicity 
clusters M3 and M13 as well as the metal poor clusters 
NGC 5466, M92 and M15. We conclude: 

• Thermohaline mixing in the presence of primordial 
enrichment can account for the carbon variation 
seen in M3 and M13. The spread of carbon along 
the RGB in both clusters can be covered by the 
same set of models. In order to match the main- 
sequence spread in M13, we require a model further 
depleted in [C/Fe]. This model however, predicts 
that some stars on the RGB should have [C/Fe] 
values lower than seen. The majority of stars in 
M13 are enriched in nitrogen whilst about half the 
stars in M3 appear to be CN strong. The results 
from the hybr id picture here are consi stent with 
the findings of iD'Antona fc Caloil (200^ and their 
requirements to match the horizontal branch mor- 
phologies of the clusters. In their work the enriched 
"^He of the second population acts as the second pa- 
rameter and helps dictate where the star falls along 
the horizontal branch. The stars in M13 that are 
most enriched in [N/Fe] suggest they were formed 
with some primary nitrogen (higher C-l-N+0 than 
those in M3). 

• Thermohaline mixing can explain carbon deple- 
tion with magnitude in NGC 5466. A single solar 
scaled model is sufficient to explain the cluster even 
though a modest CN spread has been determined 
(jShetrone et al.ll2"010D . Our models provide a good 
fit for this cluster if we adopt the LF bump deter- 
mined by Martell et al. (2008b). However, the data 
do not exclude the possibility NGC 5466 is deplet- 
ing carbon before the LF bump as seems to be the 
case for M15 and M92. 

• In M92 and M15 we have combined observations of 
[C/Fe] from multiple studies. Thermohaline mix- 
ing is unable to reproduce the evolution of carbon 
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along the giant branch in these clusters. The mod- 
els do not deplete carbon rapidly enough along the 
RGB. In addition depletion appears to begin be- 
fore first dredge- up in M92, although this may be 
due to the fact wc only have a few observations at 
low luminosity. A targeted study at low luminosity 
is required to confirm if the cluster has the same 
spread in [C/Fe] as M15. Both clusters appear to 
be mixing before the LF bump and bright stars are 
not observed near the predicted carbon abundances 
for initial [C/Fe]=0. 

• Both M92 and M15 seem to require deeper FDU in 
order for the models to fit the current observations 
but this is not the case for NGC 5466. 

• A consideration that inevitably has to serve as a 
caveat to our discussions is that the comparison 
between the thcrmohaline-mixing models and the 
[C/Fe] as a function of magnitude data may be hin- 
dered by the heterogeneity of the latter, particu- 
larly in clusters such as M92 and M15 where zero- 
point offsets of up to 0.3 dex may exist between 
the various data sources compiled. Although at- 
tempts have been made to compensate for these, 
it remains nonetheless a source of concern. Deter- 
mining the behavior of carbon to faint luminosities 
on the RGB in both M92 and M15 as weh as M3, 



using a single spectrograph and data analysis sys- 
tem could provide a more homogeneous and rigor- 
ous test of our models. Such data would need to 
include C and N values for subgiants also. In ad- 
dition, observations of lithium abundances over a 
large luminosity range can be used to better define 
the character of extra mixing, while measurements 
of nitrogen, oxygen and/or sodium abundances can 
define the patterns of primordial enrichment within 
the clusters. 
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Fig. 1. — [C/Fe] and [N/Fe] vs My for M13 (left panels) and M3 (right panels). The various studies from which the data are taken are 
listed in the Figure. In all panels the solid black curve represents our CN-weak model and the solid grey curve corresponds to our CN-strong 
model . For M13 we also provide a nitrogen enhanced model (green dashed line) and a model designed to match the main-sequence carbon 
spread of the cluster (blue dotted line). In all models the metallicity is set to Z = 0.0005 and the thermohaline mixing parameter Ct= 1000. 
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Fig. 2. — [C/Fe] vs My for NGC 5466 (top panels), M92 (middle panels) and M15 (bottom panels). In the left panels the black curve 
represents our solar scaled model (CN-weak) and the grey curve is our CN-strong model for these clusters. The solid vertical li ne repr esents 
the end of first dredge-up according to the models. The dotted line corresponds to the LF bump according to IMartell et al.' (2008 b|). Th e 
dashed lines are photometric deri vations of the LF bump for each cluster: for M92 and M15 the magnitude is taken from Nataf ct al.l 1 120111) , 
and for NGC 5466 from iFekadu et al.l I I2007I) . In the right panels wc provide lines of best fit by eye which helps highlight the strange behavior 
of the clusters. 
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Fig. 3. — f C/Fe] [N/Fe] and X(^He) vs My for M15. As above the observational data is taken from ITrefzeer et all I I1983I ') (crosses) and 
ICohen et al.l II2005I ) (diamonds). The sohd black curve corresponds to the thermohaline mixing model with standard dredge-up depth (see 
Figure |2]| where the envelope extends down to a depth of M = 0.368Mq. The grey curve represents our model where we have artificially 
homogenized beyond the depth of first dredge-up to a mass of M = O.32OM0. The soli d vertical line rep resents the end of first dredge-up 
according to the normal models. The dotted line corresponds to the LF bump of [Martell et al.l Il2008bi ) based o n a metallicity-LF b ump 
relation. This occurs at My x —0.55. The dashed line is a photometric derivation of the LF bump taken from lZoccali et al.l (|1999l ) and 
occurs at My ~ —0.42. 
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